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ABSTRACT 



Context. Change of sign of the LNRF-velocity gradient has been found for accretion discs orbiting rapidly rotating Kerr black holes 
with spin a > 0.9953 for Keplerian discs and a > 0.99979 for marginally stable thick discs. Such a "humpy" LNRF-velocity profiles 
occur just above the marginally stable circular geodesic of the black hole spacetimes. 

Aims. Aschenbach (2004) has identified the maximal rate of change of the orbital velocity within the "humpy" profile with a locally 
defined critical frequency of disc oscillations, but it has been done in a coordinate-dependent form that should be corrected. 
Methods. We define the critical "humpy" frequency v\, in general relativistic, coordinate independent form, and relate the frequency 
defined in the LNRF to the distant observers. At radius of its definition, the resulting "humpy" frequency Vh is compared to the radial 
v r and vertical v v epicyclic frequencies and the orbital frequency of the discs. We focus our attention to Keplerian thin discs and 
perfect-fluid slender tori where the approximation of oscillations with epicyclic frequencies is acceptable. 

Results. In the case of Keplerian discs, we show that the epicyclic resonance radii and r^-i (with y v : v, = 3 : 1, 4 : 1) are located 
in vicinity of the "humpy" radius r h where efficient triggering of oscillations with frequencies ~ v h could be expected. Asymptotically 
(for 1 - a < 10~ 4 ) the ratio of the epicyclic and Keplerian frequencies and the humpy frequency is nearly constant, i.e., almost inde- 
pendent of a, being for the radial epicyclic frequency v r : v h ~ 3 : 2. In the case of thick discs, the situation is more complex due to 
dependence on distribution of the specific angular momentum £ determining the disc properties. For I = const tori and 1 - a < 10~ 6 
the frequency ratios of the humpy frequency and the orbital and epicyclic frequencies are again nearly constant and independent of 
both a and I being for the radial epicyclic frequency v, : Vh close to 4. In the limiting case of very slender tori (£ ~ £ m!} ) the epicyclic 
resonance radius r 4: i ~ r h for all the relevant interval of 1 - a < 2 x 10~ 4 . 

Conclusions. The hypothetical "humpy" oscillations could be related to the QPO resonant phenomena between the epicyclic oscilla- 
tions in both the thin discs and marginally stable tori giving interesting predictions that have to be compared with QPO observations 
in nearly extreme Kerr black hole candidate systems. Generally, more than two observable oscillations are predicted. 

Key words. Black hole physics - Accretion, accretion disks - Relativity 



1. Introduction 

High frequency (kHz) twin peak quasi-periodic oscillations 
(QPOs) with frequency ratios 3:2 (and sometimes 3 :1) ar e 
observed in microquasa r s (se e, e.g.. Ivan der Klisl d2000h : 
iMcClintock & Remillardl d2004l) : iRemillardl d2005l)). I n the 
Galactic Center black hole Sgr A*. iGenzel et al.l d2003l) mea- 
sured a clear periodicity of 1020 sec in variability during a flar- 
ing event. This period is in the range of Keplerian orbital pe- 
riods at a few gravitational radii from a bl ack hole wit h mass 
M ~ 3 .6 x 10 6 M o estimated fo r Sgr A* dGhezll2004l). M ore 
recentlv lXschenbach et ail d2004l) : I Aschenbach! d2004l 120061) re- 
ported three QPO periodicities at 692 sec, 1130 sec and 2178 
sec that correspond to frequency ratios (1/692) : (1/1130) : 
(1/2178) ~ 3 : 2 : 1. Howe ver, these observational d ata are not 
quite convincing, see, e.g. dAbramowicz et al1l2004l) . In some 
galactic binary black hole and neutron-star systems, the high- 
frequency QPOs at Vhigh are accompanied with low-frequency 
QPOs at vi ow . The high-frequency and low-frequency QPOs are 
correlated and the ratio of the frequencies is observed to be 
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Vhigh ■ v\ow ~ 13 : 1 . It was first noticed by iPsaltis et al.l 
(1 1999b that the correlation between high-frequencies and low- 
frequencies exi sts f or black-hole and neu tron-star sources, later 
iMauchd d2002l) and lWarner et al.l d2003l) extended this correla- 
tion to cataclysmic variables and showed that it is obeyed by 
high-frequency quasi-coherent "dwarf nova oscillations" and the 
low-frequency "horizontal branch" oscillations. At present, there 
is no exact model explaining the ratio 13 : 1, only a qual- 
itative proposal exis ts, based on analogy wi th the 9-th wave 
from oceanography dAbramowicz et ai1l2004l) . In this concept, 
the high-frequency QPOs are connected to transient oscillatory 
phenomena at random locations in the accretion disc and are 
subject to the sid e band instability simil ar to those considered 
in oceanography dBenjamin & FeirifT967h . If a wave pulse con- 
tains initially waves of identical length and frequency Vhi g h, non- 
linearities can cause the waves with larger amplitude to move 
faster changing their wavelength. The shorter (longer) waves in 
front of (behind) the pulse cause energy to concentrate at the cen- 
ter of the pulse feeding thus the instability, the result of which 
is that every n-th wave has a higher amplitude creating low- 
frequency oscillations with frequency vi ow ~ Vhi g h/«- The value 
of n depends on details of the hydrodynamic models and it is 
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not fully understood in both o ceanography (where n ~ 9) and 
discography (where n ~ 1 2- 1 4; lAbramowicz et al. 20041). 

It was proposed bv lKluzniak & Abramow icz (2001) that the 
high frequency twin peak QPOs are r elated to the parametric o r 
forced resonance in accretion discs dLandau & Lifshitz1ll973l) . 
possibly between the radial and vertical epi cyclic oscillations 
dAliev & Gaitso^ll98ltlNowak & Lehril 19981) or the orbital and 
one of the epicyclic oscillations. Thes e oscillations could be re- 
lated to bo th the thin Keplerian discs ( lAbramowicz et al1l2003l: 
Kato 1200 lb or the thick, tor oidal accretion discs ( Rezzolla et al.l 
20031: iKluzniak et al.ll2003l) . In particular, the observations of 
high frequency twin peak QPOs with the 3 : 2 frequency ra- 
tio in microquasars can be explained by the parametric res- 
onance between the radial and vertical epicyclic oscillations, 
v v : v r ~ 3 : 2. This hypothesis, under the assumption of 
geodesic oscillations (i.e., for thin discs), puts strong limit on 
the mass-spin rela t ion for the central black hole in microquasars 
dTorok et al.ll200l lTor61dl200l: iTorok etai]l2006l) . 

lAschenbachl (120041 l200fih discovered that two changes of 
sign of the radial gradient of the Keplerian orbital veloc- 
ity as measured in the locally non-rotating frame (LNRF) 
dBardeen et al.ll 19721) occ ur in the equatorial pl ane of Kerr black 
holes with a > 0.9953. [Stuchlfk et al.l d2005l) have found that 
the gradient sign change in the LNRF-velocity profiles occurs 
also for non-geodesic motion with uniform distribution of the 
specific angular momentum ({r, 0) = const (i.e., in marginally 
stable thick discs) around extremely rapid Kerr black holes with 
a > 0.999790 The global character of the phenomenon is given 
in terms of topology changes of the von Zeipel surfaces (equiv- 
alent to equivelocity surfaces in the tori with ((r, 9) = const). 
Toroidal von Zeipel surfaces exist around the circle correspond- 
ing to the minimum of the equatorial LNRF-velocity profile, in- 
dicating possibility of development of some instabilities in that 
part of the marginally st able disc with positive gra dient of the 
orbital velocity in LNRF (IStuchlfk et al.ll2004l2005l) . 

Therefore, we consider the positive radial gradient of or- 
bital LNRF-velocity around black holes with a > 0.9953, see 
Fig. HI to be a physically interesting phenomenon, even if a direct 
mechanism relating this phenomenon to triggering the oscilla- 
tions, and subsequent linking of the oscillations to the excitation 
of radial (and vertical) epicyclic oscillations, is unknown. We 
present a basic study of the "humpy" oscillatory frequency and 
its relation to the epicyclic and Keplerian (orbital) frequencies. 
It should be stressed that recently at least two QPOs sources are 
observed, in which the rotational parameter (spin) of the central 
black hole is estimated nearly extreme, i. e., a > 0.99. Such b lack 
holes are probably observed in Sgr A* (lAschenbachll2006l) and 
in GRS 1915+105 dMcChntocketal.ll2006l) . We plan to make 
a detailed analysis of the observed frequencies and their possi- 
ble relation to the LNRF-velocity hump induced frequency and 
related epicyclic frequencies in the future work. 

lAschenbachl (12004 120061) . considering phenomena observed 
in Sgr A*, has shown that in the field of the Kerr black hole 
with a =! 0.99616, the orbit where the critical frequency sub- 
jected to the hump of the LNRF-velocity profile in such a way 
that the positive rate of change of the LNRF-velocity is max- 
imal, (v^ rit = &V (ip) /<9r) max , is located nearby r = rj,-\, where 



1 Note that the assumption of uniform distribution of the specific 
angular momentum can be relevant at least at the inner parts of the 
thick disc and that matter in the disc follows nearly geodesic circu- 
lar orbits nearby the center of the disc and in the vicinity of its inner 
edge determined by the cusp of its critical equipotential surface (see, 
lAbramowicz et alJI 19781) . 



the vertical and radial epicyclic frequencies are in the ratio of 
v v : v r = 3 : 1 and, moreover, the critical frequency v^ rit is nearly 
equal to the radial epicyclic frequency there. Undoubtly, this is 
an interesting result. However, the critical frequency introduced 
by Aschenbach is related to the rate of change of the locally mea- 
sured orbital velocity in terms of the special Boyer-Lindquist 
radial coordinate, so the coincidence v£ rit - v r obtained in this 
case is rather unrealistic. In this paper we give the critical fre- 
quency V s . , related to the maximal positive radial gradient of 
the LNRF-velocity in the "humpy" velocity profile, in the gen- 
eral relativistic, coordinate-independent form. Further, since the 
critical frequency iA t is defined locally, being connected to the 
LNRF, it has to be transformed into the form related to distant 
stationary observers, giving observationally relevant frequency 

In Section [2] we briefly summarize properties of the 
Aschenbach effect for Keplerian thin discs, and I = const 
thick discs. In Section[3] the critical frequency, connected to the 
LNRF-velocity positive gradient in the humpy profiles, is given 
in the physically relevant, coordinate independent form for the 
both Keplerian and 6 = const discs. At the radius of its defini- 
tion, the critical frequency is compared to the radial and vertical 
epicyclic frequency and the orbital frequency. In Section |4j the 
results are discussed and concluding remarks are presented. 

2. LNRF-velocity profiles of discs orbiting the Kerr 
black holes 

In the Kerr spacetimes with the rotational parameter assumed to 
be a > 0, the relevant metric coefficients in the standard Boyer- 
Lindquist coordinates read: 



gtt 



A 7 ■ 2 

A - a~ sin" i 



A sin 2 6 



gup = -■ 
S 

A' 866 



Sipip ~ £ > & 
where 

A = r 2 - 2r + a 2 , Y - r 2 + a 2 cos 2 
A = (r 2 + a 2 ) 2 - Aa 2 sin 2 6. 



2ar sin 2 6 
I ' 

= 2, 



(1) 

(2) 



(3) 
(4) 



The geometrical units, c = G = 1, together with putting the mass 
of the black hole equal to one, M = 1, are used in order to obtain 
completely dimensionless formulae hereafter. 

The locally n on-rotating frames ( LNRF) are given by the 
tetrad of 1 -forms dBardeen et al.ll 19721) 



e (D = 



/A\ 1/2 

^-j sin 6 (dip- to df), (5) 



1/2 



dr, e (fl) = S 1/2 d0, 



where 



gtip 2ar 



(6) 



(7) 



8w A 

is the angular velocity of the LNRF relative to distant observers. 
For matter with a 4-velocity U^' and angular velocity profile 
Q(r, 6) orbiting the Kerr black hole, the azimuthal component 
of its 3-velocity in the LNRF reads 



<y(<p) _ 



in 



A sin 6 

zVa 



(O - to). 



(8) 
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Fig. 1. Profiles of the equatorial orbital velocity related to LNRF in terms of the radial Boyer-Lindquist coordinate for appropriately 
chosen values of the black-hole spin a in the case of Keplerian discs (upper plots) and limiting marginally stable thick discs with 
t = Ans = const specific angular momentum distribution (lower plots). The resonant orbits r^-i and r^-\ of the epicyclic frequencies 
v v : v r together with the marginally stable (dashed line) and marginally bound (dashed-dotted line) orbits are also given. 



2.1. Keplerian thin discs 



2.2. Marginally stable tori 



In thin discs matter follows nearly circular equatorial geodetical 
orbits characterized by the Keplerian distributions of the angular 
velocity and the specific angular momentum (in the equatorial 
plane, 9 — n/2) 



£1 — £2ic(r; a) = 



= a) = 



1 



(r 3 / 2 + fl)' 

r 2 - 2ar l l 2 + a 2 
r 3/2 _ 2r i/2 + a ■ 



(9) 



(10) 



The azimuthal component of the Keplerian 3-velocity in the 
LNRF reads 



w x (r 2 + a 2 ) 2 - a 2 A - 2ar(r 3/2 + a) 



T£'(r; a) = 



r 2 (r 3 / 2 + a) VA 



(11) 



and formally diverges for r — » r+ — 1 + Vl - a 2 , where the 
black-hole event horizon is located. Its radial gradient is given 
by 

&Vf _ r 5 + a 4 (3r + 2) - 2q 3 r" 2 (3r + 1) 
dr ~ 2A 3 / 2 V7(r 3 / 2 + a) 2 

2a 2 r 2 (2r - 5) - 2ar 5/2 (5r - 9) 
2A 3 / 2 V7(r 3 / 2 + a) 2 ' 

As shown bv lAschenbachl J2004 12006). the velocity profile has 
two changes of the gradient sign (where d^V^/dr = 0) in the 
field of rapidly rotating Kerr black holes with a > a C (K) — 0.9953 
(see Fig.[T). 



(12) 



Perfect-fluid stationary and axisymmetric toroidal discs are char- 
acterized by the 4-velocity field V = (V, 0, 0, U^) with IP = 
U'(r, 6), U v = U v (r, 9), and by distribution of the specific an- 
gular momentum t = -U v /U t . The angular velocity of orbiting 
matter, Q. = U^/U', is then related to I by the formula 



Cl = — 



gtt + gup 
%tif + gipip 



(13) 



The marginally stable tori are characterized by uniform dis- 
tribution of the specific angular momentum 



t = £( r , 9) = const, 



(14) 



and are fully determined by the spacetime structure through 
equipotential surfaces of the potential W = W(r, 9) defined by 
the relations dAbramowicz et alii 19781) 



W-Wir 



In 



{V t ) h 



(U,f 



Stip 8tt8ipip 
g tt E 2 + 2g ¥ E + i 



(15) 



the subscript "in" refers to the inner edge of the disc. 
The LNRF orbital velocity of the torus is given by 



<yteO _ 



A(A - a 1 sin" 9) + \a l r L sin" 9 
I,y[A(A-2a£r)sm9 



f. 



(16) 



For marginally stable tori it is enough to consider the motion 
in the equatorial plane, 9 = n/2. Formally, this velocity vanishes 
for r — > oo and r — > r+, i.e., there must be a change of its radial 
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, inner cusp 
1 1.2 1^4 1^6 1 . 8 



Fig. 2. Von Zeipel surfaces (meridional sections). For a > fl c(T) 
and { appropriately chosen, two surfaces with a cusp, or one 
surface with both the cusps, together with closed (toroidal) sur- 
faces, exist, being located always inside the ergosphere (dashed 
surface) of a given space time. Both the outer cusp and the cen- 
tral ring of closed surfaces are located inside the toroidal equi- 
librium configurations corresponding to marginally stable thick 
discs (light-gray region; its shape is determined by the critical 
self-crossing equipotential surface of the potential W(r, 6). The 
cross (+) denotes the center of the torus. Dark region corre- 
sponds to the black hole. Figures illustrating all po ssible configu- 
rations of the von Zeipel surfaces are presented in lStuchlik et al.l 
(12005b . Here we present the figure plotted for the parameters 
a = 0.99998, I = 2.0065. Critical value of the von Zeipel radius 
corresponding to the inner and the outer self-crossing surface is 
^c(in) — 3.429 and 7? C (out) - 3.804, respectively, the central ring 
of toroidal surfaces corresponds to the value Reenter - 3.817. 
Interesting region containing both the cusps and the toroidal sur- 
faces is plotted in detail at the left lower figure. Right lower fig- 
ure shows the behaviour of the von Zeipel radius in the equatorial 
plane. 



gradient for any values of the parameters a and {, contrary to 
the case of Keplerian discs. The radial gradient of the equatorial 
LNRF velocity of I = const tori reads 

&Vf [[A + (r- l)r][r(r 2 + a 2 ) - 2a{f - a)] 



[r(r 2 + a 2 ) - 2a(t - a)] 2 VA 
r(3r 2 + a 2 )A } 
[r(r 2 + a 2 ) - 2a(( - a)] 2 VA 



(17) 



so it changes its orientation at radii determined for a given I by 
the condition 



= 4x(r; a) = a + 



r 2 [(r 2 + a 2 )(r- 1) - 2rA] 
2a[A + r{r - 1)] ' 



(18) 



Of course, for both thick tori and Keplerian discs we must 
consider the limit on the disc extension given by the innermost 
stable orbit. For Keplerian discs this is the marginally stable 
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Fig. 3. Velocity difference AT** = < V^ X - f V^ 1 as a function 
of the rotational parameter a of the black hole for both the thin 
(Keplerian) disc and the marginally stable (non-Keplerian) disc 
with I = 4, s . 



geodetical orbit, r m « r ms , while for thick tori this is an un- 
stable circular geodesic kept stable by pressure gradients and lo- 
cated between the marginally bound and the marginally stable 
geodetical orbits, r m b < r; n < r ms , with the radius being deter- 
mined by the specific angular momentum t = const e (Z ms , Z m b) 
through the equation I = €jt(r\ a); £ ms (£ m b) denotes specific an- 
gular momentum of the circular marginally stable (marginally 
bound) geodesic. 

Detailed discussion o f IStuchKk et all (l2005h shows that two 
physically relevant changes of sign of jdr in the tori oc- 
cur for Kerr black holes with the rotational parameter a > 
fl c(T) 0.99979 (see Fig. [TJ- The interval of relevant val- 
ues of the specific angular momentum I € {£ ms (a), { ex (max)(a)), 
where 4x(max)(«) corresponds to the local maximum of the func- 
tion (1181) . grows with a growing up to the critical value of 
«c(mb) — 0.99998. For a > a c ( m t>), the interval of relevant values 
of I 6 ({ ms (a), ( m b(a)) is narrowing with the rotational parame- 
ter growing up to a = 1, which corresponds to a singular case 
where i ms {a — 1) — { m b(a — 1) = 2. Notice that the situation 
becomes to be singular o nly in terms of the sp ecific angular mo- 
mentum; it is shown (see lBardeen et al.l 19721) that for a — 1 both 
the total energy E and the axial angular momentum L differ at 
r ms and r m b, respectively, but their combination, I = L/E, giving 
the specific angular momentum, coincides at these radii. 

It should be stressed that in the Kerr spacetimes with a > 
a C (T), the "humpy" profile of < V)p(r; a) occurs closely above the 
center of relevant toroidal discs, at radii corresponding to stable 
circular geodesies of the spacetime, where the radial and vertical 
epicyclic frequencies are also well defined. 

A physically reasonable way of defining a global quantity 
characterizing rotating fluid configurations in terms of the LNRF 
orbital velocity is to introduce, so-called, von Zeipel radius de- 
fined by the relation 



n = 



dp) 



= (1 -cj£)q, 



(19) 



which generalizes in a nother way as compared with 
dAbramowicz et all 1 1995b the Schwarzschildian definition 
of the gyration radius g dAbramowicz et all 119931) . Note that, 
except for the Schwarzschild case a = 0, the von Zeipel 
surfaces, defined as the surfaces of 7?(r, 9; a, l) = const, do not 
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coincide with those introduced by iKozlowski et al.l d!978l) as 
the surfaces of constant ^/Q.Q 

In the case of marginally stable tori the von Zeipel sur- 
faces % = const coincide with the equivelocity surfaces 
T^Cr, 0; a, t) = < V% ) = const. Topology of the von Zeipel 
surfaces can be directly determined by the behaviour of the von 
Zeipel radius in the equatorial plane 



R(r, = tt/2; a, €) 



r(r 2 + a 2 ) - 2a(( - a) 
rVA 



(20) 



The local minima of the function d20b determine loci of the cusps 
of the von Zeipel surfaces, while its local maximum (if it ex- 
ists) determines a circle around which closed toroidally shaped 
von Zeipel surfaces are concentrated (see Fig.[2]i. Notice that the 
minima (maximum) of !R(r, 9 = n/2; a, C) correspond(s) to the 
maxima (minimum) of "Vljf \r, 6 - n/2; a, €), therefore, the in- 
ner cusp is always physically irrelevant being located outside of 
the toroidal configuration of perfect fluid. Behaviour of the von 
Zeipel surfaces nearby the center and the inner edge of the thick 
discs orbiting Kerr black holes with a > a c <j) = 0.99979, i.e., 
the existence of the von Zeipel surface with an outer cusp or the 
surfaces with toroidal topology, suggests possible generation of 
instabilities in both the vertical and radial direction. 



2.3. Velocity profiles with a hump 

Behavior of <v!jf ) (r; fl > i) and Vf(n a) is illustrated in Fig. [TJ 
With a growing in the region of a e (a C (T), 1) (a e (a C (K> 1)), 
the difference A^V ( f = <V W 



T(max) 



-<V, 



if) 

T(min) 



^ K _ K K(max) 



"Vwrojjjp grows (Fig. [3} as well as the difference of radii, Arj = 
(Ar K = ?"K(max) - '"K(min)), where the local extrema 

of <y^) (VyM) occurj see p ig g] 

In terms of the redefined rotational parameter (1 -a), the 
"humpy" profile of the LNRF orbital velocity of marginally sta- 
ble thick discs occurs for discs orbiting Kerr black holes with 
(1 - a) < 1 - a C (T) — 2.1 X 10~ 4 , which is more than one order 
lower than the value 1 - a C (K) — 4.7 X 10~ 3 found bv lAschenbachl 
d2004l) for the Keplerian thin discs. Moreover, in the thick discs, 
the velocity difference A'Vjf is smaller but comparable with 
those in the thin discs (see Fig. [3]). In fact, we can see that for 

(T) 

jg - 0.07. 



a 



1, the velocity difference in the thick discs A*y^ 
while for the Keplerian discs it goes even up to A*VjS 



0.02, 



3. Humpy frequency and its relation to epicyclic 
frequencies 

In Kerr spacetimes, the frequencies of the radial and latitudinal 
(vertical) epicyclic oscillations related to an equatorial Keplerian 
circular orbit at a given r are determined by the formulae (e.g., 
lAliev & Gal tsov 198lt lNowak & Lehril 19981) 



v 2 = y|(l - 6r l + 8ar V2 - 3a 2 r 2 ), 
v 2 v = v 2 = 4(1 - 4ar' /2 ■ 



+ 3a 2 r 2 ), 



(21) 
(22) 



where vk = Qk/2tt. A detailed analysis of properties of 
the epi cyclic frequencies can be found in iTorok & Stuchlfk! 
d2005alfbl) . The epicyclic oscillations with the frequencies 
v r , v v can be related to both the thin Keplerian discs 
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a 

(b) 

Fig. 4. Positions of local extrema of r V (jp ) ( m g_L coordinates) 
for Keplerian discs (a) and marginally stable discs with ( — £ ms 
(b) together with the locations of resonant orbits ry\ and r^ : \ 
(where the resonance between the vertical and radial epicyclic 
oscillations takes place) in dependence on the rotational param- 
eter a of the black hole. Dashed curve corresponds to the max- 
imum positive values of the LNRF orbital velocity gradient in 
terms of the proper radial distance where the critical frequency 
v^ rit is defined, boundaries of shaded regions correspond to orbits 
where the velocity gradient giving the characteristic frequency, 
tfV^ldR, reaches (a) 99%, 90%, 80% and (b) 99%, 95%, 90% 
of its maximum. 



and thick, toroidal 



dAbramowicz & Kluzniafll2000l: lKatoll2005h and 
discs dRezzolla etaL 2003; Kl uzniak et all2004h . 

lAschenbachl d2004l 120061) defined the characteristic (criti- 
cal) frequency of any related mechanism possibly exciting the 
disc oscillations in the region of positive gradient of its LNRF- 
velocity < V^- ) by the maximum positive slope of "V^: 



v A - 



Br 



(23) 



- For more details see IStuchliketaLl f2005). 



This frequency has to be determined numerically and we have 
done it for both the Keplerian discs and the marginally stable 
discs with £ = £ ms = const, see Fig.[5]and TableQ] 

Although there is no detailed idea on the mechanism gen- 
erating the "hump-induced" oscillations, it is clear that the 
Aschenbach proposal of defining the characteristic frequency 
deserves attention. It should be stressed, however, that a detailed 
analysis of the instability could reveal a difference between the 
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Fig. 6. Determination of the critical "humpy" frequency, (a) Positive parts of the "coordinate" and "proper" radial gradient d'V^ /dr 
and d^V^'/dR for a given value of the rotational parameter a in the Keplerian disc, (b) Proper radial distance of the loci of 
{d'V^ldR)^ measured from the marginally bound orbit for both the Keplerian discs (Rf rit ) and I — { ms perfect-fluid tori (/? c ™ s t ). 
Proper radial distance to the marginally stable orbit (R m $) is also shown. 



characteristic frequency and the actual observable one, as the lat- 
ter should be associated with the fastest growing unstable modfl 
Moreover the frequency vjL , defined by Eq. d23l ). represents an 
upper limit on the frequencies of the hump-induced oscillations, 
as it is given by maximum of the LNRF-velocity gradient in the 
humpy part of the velocity profile. 

In the following we assume that the characteristic (critical) 
frequency is a typical frequency of oscillations induced by the 
conjectured "humpy instability", and that the humpy oscillations 
could excite oscillations with the epicyclic frequencies or some 
combinational frequencies, if appropriate conditions for a forced 
resonance are satisfied in vicinity of the radius where the humpy 
oscillations occur. 

In situations where the general relativity is crucial, it is nec- 
essary to consider d^V^'/dR, where R is the physically relevant 
(coordinate-independent) proper radial distance, as this is an ap- 
propriate way for estimating the characteristic frequencies re- 
lated to local physics in the disc. Then correct general relativis- 



tic definition of the critical frequency for possible excitation of 
oscillations in the disc is given by the relations 



*cnt 



dR 



dR 



grrdr 



(24) 



where <V ( ^ = < V^ ) (r; a) in thin Keplerian discs, and "V** = 

f Vj(r, I, a) in marginally stable thick discs. Of course, such a 
locally defined frequency, confined naturally to the observers or- 
biting the black hole with the LNRF, should be further related 
to distant stationary observers by the formula (taken at the B-L 
coordinate r corresponding to (cKV^' / 'dR) mdx ) 



(25) 



We thank to the referee for pointing out this possibility. 



We suggest to call such a coordinate-independent and, in princi- 
ple, observable frequency the "humpy frequency", as it is related 
to the humpy profile of *V^, and denote it Vh- Again, the physi- 
cally relevant humpy frequency Vh = v£,, connected to observa- 
tions by distant observers and exactly defined by Eqs. ( |24T i and 
(l25"l l. represents an upper limit on characteristic frequencies of 
oscillations induced by the hump of the LNRF-velocity profile, 
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Fig. 7. Spin dependence of the ratio of the radial epicyclic fre- 
quency and the "humpy frequency" related to distant observers. 
The ratio is given in the radius of definition of the humpy fre- 
quency rh. In the interval of 1 - a e (1.7 x 10~ 3 , 10~ 4 ), the ratio 
rapidly falls down, to the asymptotic value of 3 : 2 starting at 
a ~ 10~ 4 . Then an exact 1/M scaling holds with frequencies 
depicted in the figure. Notice that at the Aschenbach's value of 
a =i 0.99616, for which the resonant orbit with v v : v r ~ 3 : 1 
is close to rh, there is v v /vh ~ 12, analogous to the ratio of high 
and low frequency QPOs. 



and the realistic humpy frequencies, as observed by distant ob- 
servers, can be expected close to but smaller than vfj,. Further, we 
denote the B-L radius of definition of the humpy oscillations 
frequency, where dV^/dR = (d<V (ip) /dR) m - dx . Of course, in re- 
alistic situations the hump-induced oscillation mechanism could 
work at the vicinity of rj,, with slightly different frequencies; we 
should take into account that the shift of the radius, where the 
mechanism works, shifts both the locally measured (LNRF) fre- 
quency (Eq. d24t ) and the frequency related to distant observers 
(Eq. (T2"5ll). The zones of radii, where the critical frequency v^ rf 
differs up to 1%, 10% and 20% of its maximal value (given by 
(ST^/d^W) for thin (Keplerian) discs or 1%, 5% and 10% of 
its maximum for marginally stable discs with I = { ms , are given 
in Fig. |U 

An analogical relation to Eq. dZ5b can be written also for 
the Aschenbach critical frequency iA ( , giving the Aschenbach 
frequency related to distant observers vQ,. Because the velocity 
gradient related to the proper distance R is suppressed in com- 
parison with that related to the Boyer-Lindquist coordinate dis- 
tance r, there is y* , < v\. The situation is illustrated in Fig.B] 

' r-P — 

Moreover, Fig. [6] shows mutual behaviour of the "coordinate" 
and "proper" radial gradient d'V^ jdr and <9<yw jdR in the re- 
gion between the local minimum and the outer local maximum 
of the orbital velocity < V^ of t — £ ms = const discs for an ap- 
propriately chosen value of the rotational parameter a. It is inter- 
esting to compare the Aschenbach frequencies (defined in terms 
of the B-L coordinate r) with the critical frequencies defined in 
terms of the proper radial distance R. Characteristic frequencies 

v crif v °°> ^oit' ^ gi ven m Table[T]for some typical values of 
the rotational parameter a for both Keplerian discs and limiting 
t = const tori with £ = i ms . 

The physically and observationally relevant frequency con- 
nected to the LNRF- velocity gradient sign change is given by the 
frequency Vh = v£, corresponding to the locally "hump-induced" 
oscillations taken from the point of view of distant stationary 
observers. In order to obtain an intuitive insight into a possible 
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Fig. 8. Spin dependence of the ratios of the radial (v r ) and verti- 
cal (v v ) epicyclic frequencies, and the Keplerian frequency (vk) 
to the thin-disc humpy frequency related to distant observers 
(vh). Further the ratio of the epicyclic frequencies is given at the 
radius of definition of the humpy frequency. All the frequency 
ratios are asymptotically (for 1 - a < 10~ 4 ) constant. There si 
vk : v v : v r : Vh ~ 46 : 1 1 : 3 : 2. Therefore, we can expect 
some resonant phenomena on the ratio of v r : Vh ~ 3 : 2, and 



4 that could be both correlated. 




log(l - a) 

Fig. 9. Interval of humpy frequencies for the marginally stable 
thick discs with I e (£ ms , ^ m b) as a function of the black-hole 
spin a. For a — > 1, the interval is narrowing and asymptoti- 
cally reaching the value of 150Hz(M/M o ) _I . Dotted curve cor- 
responds to the humpy frequencies of marginally stable slender 
tori with I = £ C rit, for which the critical von Zeipel surface con- 
tains two cusps (as it is demonstrated for one special case in the 
left panel of the figure; the thick part of the torus is given by the 
light-gray region). 



observational relevance of Vh, it is useful to compare it with the 
frequencies of the radial and vertical epicyclic oscillations, v r 
and v v , and the orbital frequency of the disc, v or b = Q./2n, where 
Q. is given for both thin and thick discs by Eq.[l3]and the appro- 
priate distribution of the specific angular momentum 6. The most 
interesting and crucial phenomenon is the spin independence of 
the frequency ratios for extremely rapid Kerr black holes. The 
results are given in Figs. [711101 Further we can see (Figs. |4j> that 
the resonant epicyclic frequencies radii r^i and are located 
within the zone of the hump-induced oscillation mechanism in 
both thin discs and marginally stable tori. 
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We would like to call attention to the fact that in Keplerian 
discs the sign changes of the radial gradient of the orbital veloc- 
ity in LNRF occur nearby the r — ry \ orbit (with v v : v r = 3 : 1), 
while in the vicinity of the r = ryz orbit (with v v : v, • = 3 : 2), 
d^V^/dr < for all values of a for both Keplerian discs and 
marginally stable tori with all allowed values of £. The para- 
metric resonance, which is the strongest one for the ratio of the 
epicyclic frequencies v v : v r = 3 : 2, can occur at the r = r^a 
orbit, whil e its effect is much smalle r at the radius r = r^\, as 
noticed by lAbramowicz et al] d2003l) . Nevertheless, the forced 
resonance may take place at the r^i orbit. Notice that the forced 
resonance at r = ry_] can generally result in observed QPOs fre- 
quencies with 3:2 ratio due to th e beat frequencies allowe d for 
the forced resonance as shown in lAbramowicz et al I (12004 . But 
the forced resonance at ry t \ between the epicyclic frequencies, 
induced by the humpy profile of ( V (fp \ seems to be irrelevant 
in the case of microquasars, since all observed frequencies lead 
to the values o f the r otational parameter a < a c (K), as shown 
by iTorok et al] (120051) . On the other hand, the LNRF- velocity 
hump could induce the forced resonance between another (non- 
epicyclic) frequencies as well, and thus being relevant also for 
microquasars lik e the nearly extreme Ke rr black hole candidate 
GRS 1915+105 dMcClintock et alj|2006l) . 

The marginally stable tori have a structure that depends on 
the value of the specific angular momentum t e (( ms , ( m b)- The 
oscillations of slender tori {I » £ ms ) have frequencies equal to 
the epicyclic frequencies relevant for test particle motion, but 
the frequencies of non-slender tori are different, as shown for 
pseudo-Newtonian tori (Sramkova 2 0051 IMaes et al.ll2006l) and 
expected for tori in the strong gravitational field of Kerr black 
holes. Therefore, comparison of the humpy frequencies and the 
epicyclic frequencies is relevant for the slender tori only. 

The humpy frequency is defined for all a > 0.99979 
and all £ e ({ ms , £ m b), see Fig. [9] It is important that in 
the field of Kerr black holes with 1 - a < 10 8 , there is 
Vh(a, () 150 Hz (M/M ) _1 independently of a and t. Further, 
it is shown that physically important case of tori admitting evo- 
lution of toroidal von Zeipel surfaces with the critical surface 
self-crossing in both the inner and the outer cusps is allowed at 
£ — 4rft> where £ c rit > £ ms only slightly differs from ^ ms , i.e., 
such tori can be slender, see Fig. [9] The ratios of v r /vh, v v /vh 
and v /vh are given for the tori with 6 ~ £ ms in Fig. \W\ Their 
asymptotical values, valid for 1 - a < 10~ 6 , are independent of 
both a and £. 



4. Concluding remarks 

The equality of iA. and v r for the Kerr black holes with a 
0.99616, indicating direct relation of the Aschenba ch character- 
istic frequen cy and the radial epicyclic frequency dAschenbachl 
12004 120061) . is rather only an accidental coincidence, because 

is defined in a coordinate-dependent way. The physically rel- 
evant frequency Vh = vf, cannot be directly related to the radial 
epicyclic frequency in Keplerian discs, as < v r for all relevant 
values of a e (a C (K), 1)- Nevertheless, the behaviour of the ratio 
v r /vf, indicates some interesting consequences (see Hig. [7J. 

First, for thin (Keplerian) discs around the Kerr black holes 
with a 0.99616, when the ratio of epicyclic frequencies v v : 
v r ~ 3 : 1 at the radius of definition of vf,, we find v r : ~ (11- 
13) : 1, i.e., in such a situation the frequency induced by the 
positive gradient of the LNRF-velocity profile could be related 
to the low-frequency oscillations. However, such explanation is 
restricted to extremely rapidly rotating black holes and, contrary 
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Fig. 10. Spin dependence of the ratios of the radial (v r ) and ver- 
tical (v v ) epicyclic frequencies, and the orbital frequency (v ) of 
the marginally stable I - £ ms disc to the thick-disc humpy fre- 
quency related to distant observers (vh). All the frequency ratios 
are asymptotically (for 1 - a < 10~ 6 ) almost constant. 



to the idea of the 13-th wave (lAbramowicz et al" 1 12004b . cannot 
be extended to other black-hole, neutron-star, and white-dwarf 
systems. Therefore, this has to be taken as a kind of curiosity 
working for a very special class of black-hole systems only. 

Second, for thin (Keplerian) discs around the Kerr black 
holes with a > 0.9999, there is the ratio of v r : v£, ~ 3 : 2, 
and v v : vf, ~ 11 : 2, independently of a. Assuming that 
the oscillations at the humpy frequency Vh = v£, could be re- 
ally directly detected by distant observers, for such black holes 
with 1 - a < 10 4 the high-frequency twin peak QPOs with 
3 : 2 ratio could be explained independently of the standard res- 
onant phenomena, if we focus on the asymptotic behaviour of 
v r : v^, ~ 3 : 2. Moreover, for such extremely rapid Kerr black 
holes with 1 - a < 10~ 4 , we could consider triples of frequencies 
taken in rational ratios v v : v, ■ : v£, ~ 1 1 : 3 : 2, if the epicyclic 
oscillations are excited by the LNRF-velocity hump. Such fre- 
quency ratios could be observed mainly in disc systems around 
supermassive black holes in galactic nuclei that are expected to 
be extremely fast rotating; especially Sgr A* should be tested 
very carefully for this possibility. For Kerr black holes with the 
spin parameter 1 - a > 10~ 4 , the frequency ratio is different and 
depends strongly on the spin a (see Figs. [7] and [8]). 

Considering also the Keplerian frequency we find the ratio 
of vk : having a local minimum for a 0.99965 and a nearly 
constant value for 1 - a < 10~ 5 , where vk : ~ 23. In the 
field of Kerr black holes with 1 - a < 10~ 5 , the frequency ratios 
vk : v v : v r : vf, ~ 46 : 1 1 : 3 : 2 are almost independent of a. 
Thus for the extremely rapid Kerr black holes the 1 jM scaling 
of considered frequencies is quite exact. Note that in such a case 
there is v r : Vh ~ 3 : 2 and the ratio vk : v v is close to the ratio 4: 1 
at the radius of definition of the humpy frequency. This indicates 
a possibility of "doubled" resonant phenomena with the special 
frequency ratios in Keplerian discs orbiting extremely rapid Kerr 
black holes (1 -a < 10~ 5 ). 

Third, the hump-induced oscillations with frequencies Vh < 
vf, could be generated in a zone around h, (vh — v^, at h,), 
where the resonant phenomena between the radial and vertical 
epicyclic oscillations could enter the game, namely at the ratios 
of v v : v r = 3 : 1 and 4:1. Interesting resonant phenomena 
could be then expected when the v r : Vh corresponds to the ratio 
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Table 1. Characteristic frequencies in units of (M/M Q ) 1 Hz (M/M Q is the mass of the Kerr black hole in units of mass of the 
Sun), corresponding to critical frequencies iA ( , v£,, v^ rit , vf, defined in the text, are given for appropriate values of the black-hole 
spin. Maximal values of the frequencies related to the stationary observer at infinity are bold-faced. Note that only the frequencies 
vf, have physical meaning for direct comparison with the frequencies of orbital oscillations v or b, Vr, v v related to the observers at 
infinity. 
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of small integer numbers. Especially the case of v r : Vh ~ 3 : 2 
in spacetimes with 1 - a < 10~ 4 is worth of attention. In general, 
observationally relevant should be the resonances represented 
b y frequency ratios in smal l integer numbers p : q. As shown 
in lLandau & Lifshitzl (1 19731) . the relevance of resonant phenom- 
ena depends on the order of resonance n = max(p, q), and falls 
steeply (in powers) with increasing value of n; in fact they argue 
that relevant resonant phenomena could be expected for n < 4. 
Therefore, the frequency ratios such as 23:1, 11:2, 11:3 appear 
to be quite irrelevant in realistic resonance models. 

Recall that there is a well known Thorne limit giving the 
maximum spin of the Kerr black hole in systems with thin ac- 
cretion discs, a max a 0.998, determined by the back-reaction of 
photons radiated from the disc and c aptured by the black hole 
dPage & Thornelll974l: iThornd [19741) . If the hump-induced os- 
cillations and related epicyclic frequencies will be observed in 
ratios corresponding to the asymptotic region of a > 0.9999 
for Keplerian discs, the Thorne model should be corrected, 
e.g., by effect of an occultation of the disc. In the case for 
which the Thorne limit turns out to be realistic, the hump- 
induced oscillations have to be restricted on the spin interval 
a € (0.9953, 0.998). We expect the Thorne limit being relevant 
for smooth thin discs, while the overcoming of a max could be 
expected in highly turbulent discs with toroidal internal parts. 

For thick discs the situation is much more complex, being de- 
pendent on both the rotational parameter (spin) a and the specific 
angular momentum £. The range of maximal humpy frequencies 
for a given spin a is plotted in Fig. [9] and is determined by their 
evaluation in limiting values of the specific angular momentum 
I relevant for the "humpy" effect in marginally stable thick ac- 
cretion discs (see the discussion in Sec. 12. 2b . The minimal value 
corresponds to £ ms {a) while the maximal value, in dependence 
of a, corresponds to 4 x(max) (a) (for 0.99979 < a < 0.99998) 
or £ mh (a) (for 0.99998 < a < 1). Notice that asymptotically 
(for 1 - a < 10~ 8 ) both Vh(ms) an d Vh(mb) coincide on the line of 
150 Hz (M/M ) _1 . Clearly, the same is true for the humpy fre- 
quencies related to discs with any relevant £ e {£ ms , £ m b)- The 
spin-dependence of the ratio of the humpy frequency and the 
epicyclic and orbital frequencies (taken at the radius of definition 
of the humpy frequency) for the case of limiting £ = £ ms discs is 
given in Fig. [10] Again we obtain asymptotically constant (spin- 
independent) ratios for black holes with 1 - a < 10~ 6 , where 
v r : v h(ms) ~ 15 : 4, v v : v r = 4.39, v v : v h(ms) = 16.54, v orb : 



v, = 25.4, and v or b : Vh — 95.4. It should be stressed that for 
the holes with 1 - a < 10~ 6 the same ratios with the humpy fre- 
quencies are obtained for the discs with any £ e (£ ms , £ m b), as 
£mb — * Am — * 2 for a — ¥ 1 . The asymptotically constant values 
of the frequency ratios correspond to the rational value only in 
the case of v, : vj, ~ 15 : 4. Of course, we could find some ra- 
tional ratios for any ones, if 1 - a > 10~ 6 . On the other hand, we 
directly see (Fig.[4j5) that for the very slender marginally stable 
tori (£ as £ ms ) the resonant phenomena on epicyclic frequencies 
with v v : v r = 4 : 1 ratio appear in very close vicinity of the 
humpy radius r^, making thus a very special prediction on the 
QPOs frequencies observed in such hypothetical systems with 
Kerr black holes having 1 - a < 2 x 10~ 4 . In the marginally 
stable slender tori the resonant phenomena between the radial 
epicyclic and humpy oscillations, taking place at the humpy ra- 
dius rh, and between the vertical and radial epicyclic oscillations 
near the humpy radius, both with the ratio ~ 4 : 1, could be ob- 
servationally relevant only, but their relevance is expected to be 
lower than that of the frequency ratios 3:2 and 3:1 in Keplerian 
discs. 

Finally, it should be stressed that at present no direct mech- 
anism triggering the LNRF velocity hump excited oscillations is 
known, being a challenge for investigation, since the existence 
of the toroidal von Zeipel surfaces (see Figs.[2][9]l brings some 
indication of possible triggering of instabilities in both radial and 
vertical directions leading to oscillations in accretion discs. The 
predictions for the ratio of the humpy and epicyclic or Keplerian 
(orbital) frequencies presented here for both thin discs and slen- 
der tori have to be compared with observations made in nearly 
extreme Kerr black hole systems. In the case of the humpy oscil- 
lations excited systems we could observe more than two QPOs 
with frequencies in rational ratio. It seems that in the X-ray vari- 
able binary system (microquasar) with the nearly extreme Kerr 
black hole candidate GRS 1915+105 four oscillations with re- 
lated frequencies have been observed, what brings a large field 
for testing the predictions of the "LNRF-velocity hump excited 
oscillations" model. The tests have to be done in a close con- 
nection to both the related resonance model and the results of 
the spectral ana lysis of the X-ray continuum, as obser ved in 
GRS 1915+105 jRemillardll2005HMcClintock et alj2006l) Fl We 



4 If the results of I Aschenbach et all d2004h will be confirmed by 
more precise observations, the Galactic Center black hole system 
Sgr A* could serve as another example of the nearly extreme Kerr black 
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believe that a synergy effect of such studies could lead to deeper 
understanding of X-ray binary systems, namely microquasars. 
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